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Abstract 

We report an improved measurement of the D'^D^ and D^D^ decays based 

on 656.7 x 10^ BB events collected with the Belle detector at KEKB. We measure the branching 

fraction and charge asymmetry for the B+ -> D+D^ decay: B{B+ D+D^) = (3.85±0.31±0.38) X 
10~^ and Acp{B^ D^D^) = 0.00ib0.08±0.02, where the first error is statistical and the second is 
systematical. We also set the upper limit for the 5° D^D^ decay: B{B^ -> D^D^) < 0.42 x 10"^ 
at 90% CL. 

PACS numbers: 13.25.Hw, 14.40.Lb 
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INTRODUCTION 



Recently, evidence of direct CP violation in —>■ D^D^ decays was observed by the 
Belle collaboration [l|. A possible interpretation is that there is a sizable penguin contribu- 
tion to this decay. If this is the case, a similar effect should be seen in the charged mode 

D^D^, and has been already observed by Belle 0] and confirmed by BaBar 
The inclusion of charge conjugate states is implicit throughout this paper. In this paper, 
we report an improved measurement of the branching fraction and charge asymmetry for 
— s> D^D^ decay and we also search for the B^ D^D^ decay. The latter can only be 
produced by W exchange diagram. We use a data sample of 656.7 x 10^ BB events collected 
with the Belle detector at the KEKB collider [i]. 

BELLE DETECTOR 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [Ef. For the first sample of 152 million 
BB pairs, a 2.0 cm radius beam pipe and a 3-layer silicon vertex detector were used; for 
the latter 505 million BB pairs, a 1.5 cm radius beam pipe, a 4-layer silicon detector and a 
small-cell inner drift chamber were used 

EVENTS SELECTION 

The track transverse momentum is required to be higher than 0.075 GeV/c in order to 
reduce combinatorial background. 

For charged particle identification (PID), the measurement of the specific ionization 
(dE/dx) in the CDC, and signals from the TOF and by ACC are used. Charged kaons 
are selected with PID criteria that have an efficiency of 88% with a pion misidentification 
probability of 8%. All charged tracks that are consistent with a pion hypothesis that are 
not positively identified as electrons are treated as pion candidates. 

Neutral kaons are reconstructed in the decay Ks — > h'^tt'; no PID requirements are ap- 
plied for the daughter pions. The two-pion invariant mass is required to be within 9 MeV/c^ 
(~ 3cr) of the mass and the displacement of the tt+tt^ vertex from the IP in the trans- 
verse (r — (p) plane is required to be between 0.2 cm and 20 cm. The Ks momentum and 
the vector from the IP to the n'^ir' vertex are required to be collinear in the r — ip plane to 
within 0.2 radians. 

Photon candidates are selected from ECL showers not associated with charged tracks. 
An energy deposition of at least 75 MeV and a photon-like shape of the shower are required 
for each candidate. A pair of photons with an invariant mass within 12 MeV/c^ (~ 2.5a) of 
the 71^ mass is considered as a 7r° candidate. 

We reconstruct mesons in the K^tt^, K^-k^-k^'k^ and K^ti^-k^ decay channels. The 
candidates are reconstructed in K^-k^-k^ and Ks^^^ final states. We require the invariant 
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mass of the D^{D^) candidates to be within 11 MeV/c^ {1.5a for K~tt^tt'^ and 2.5a for 
other modes) of the nominal D^{D~^) mass. For — K'lr'^n'^ decay, we require that 
the 71^ momentum be greater than 0.35 GeV/c in order to reduce further the combinatorial 
background. To suppress the large background from D^D^ decays with the 

from misreconstructed as a pion, none of the pions from should be consistent with 
the kaon hypothesis. This requirement has an efficiency of 93% and kaon misidentification 
probability of 9%. 

We combine and {D^) candidates to form (B^) candidates. These are identi- 
fied by their center-of-mass (CM) energy difference, AE = Ei) — -Ebeam, and the beam 
constrained mass, Mbc = ^/ E^eam ~ ^K^iPiYi where -Ebeam is the beam energy and Pi and 
Ei are the momenta and energies of the decay products of the B meson in the CM frame. 
We select events with Mbc > 5.2 GeV/c^ and \A.E\ < 0.3 GeV, and define a B signal region 
of \AE\ < 0.02 GeV, 5.273 GeV/c^ < Mbc < 5.287 GeV/c^. In an event with more than one 
B candidate, we choose the one with smallest from D mass fit. We use a Monte Carlo 
(MC) simulation to model the response of the detector and determine the efficiency 

Variables that characterize the event topology are used to suppress background from the 
jet-like e+e~ — >• qq continuum process. We require | cos6'thr| < 0.8, where 6'thr is the angle 
between the thrust axis of the B candidate and that of the rest of the event; this condition 
rejects 77% of the continuum background while retaining 78% of the signal. To suppress 
high background in the D^D^ final state, we use a Fisher discriminant, JF, that is based 
on the production angle of the B candidate, the angle of the B candidate thrust axis with 
respect to the beam axis, and nine parameters that characterize the momentum fiow in the 
event relative to the B candidate thrust axis in the CM frame [8|. We impose a requirement 
on JF that rejects 52% of the remaining continuum background and retains 86% of the signal. 



RESULTS 

Measurement of the branching fractions 

The AE and Mbc distributions for 5+ D+D^ and 5° ^ D^D° candidates are pre- 
sented in Figs. [1] and [2l The region AE < —0.1 GeV is excluded from the fit to avoid 
contributions from B —>■ DD* decays. 

The signal probability density function (PDF) is described by double Gaussian for AE 
and a single Gaussian for M^c. The Aii^-Mbc correlation is taken into account. We use 
the 5+ D+D^ events in our data sample to calibrate the means and resolutions of the 
signal shape. The continuum, BB and B DgD^ background contributions are described 
separately. We use a linear function for AE and threshold function for Mbc to describe 
the continuum PDF. The BB background is modeled by a quadratic polynomial for AE, 
a threshold function for Mbc combined with a small peaking component (wide Gaussian on 
AE, Gaussian for Mbc). The shape of peaking background and threshold function param- 
eters are fixed from the generic BB MC. The AE linear slope and quadratic term are free 
parameters. The peak in the AE distribution near —70 MeV coming from the B DgD^ 
decay is described by a Gaussian for AE and a Gaussian for Mbc. Again, we use B^ D^D^ 
to obtain the parameters of this PDF. 

We determine the signal yield from the two-dimensional (2D) fit to the Aii^-Mbc distri- 
bution. As a cross-check, we also do separate one-dimensional fits for the AE and Mbc 
distributions, with the value of the other variable being in the signal region. The results are 
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FIG. 1: AE (a) and M^c (b) distributions for the D^D^ candidates. Each distribution 

is the projection of the signal region of the other parameter. Points with errors represent the 
experimental data, crosshatched histograms show the BB MC and curves are projections from the 
two dimensional fits. 




FIG. 2: AE (a) and Mbc (b) distributions for the 5° D^D^ candidates. Each distribution 
is the projection of the signal region of the other parameter. Points with errors represent the 
experimental data, crosshatched histograms show the BB MC and curves are projections from the 
two dimensional fits. 



TABLE I: Yields from the AE, Mbc and 2D (AE'-Mbc) fits, detection efficiencies, and corresponding 
branching fractions. Errors are statistical only. Upper limits are at the 90% confidence level. 



Decay channel 


AE yield 


Mbc yield 


2D yield 


e, 10-4 


B, 10-4 




-4.5 ±29.7 


5.7 ± 28.6 


0.4 ±24.8 (< 41) 


16.4 


< 0.38 




366.4 ± 31.8 


376.4 ± 30.7 


369.7 ± 29.4 


14.6 


3.85 ±0.31 
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TABLE II: Charged B mesons yields from the AE, Mbc and 2D (A£J-Mbc) fits. Errors are 
statistical only. 



Decay channel 


AE yield 


Mbc yield 


2D yield 


B+ D+D° 


183.9 ±21.5 


184.4 ±21.4 


184.2 ± 20.4 


B- D-D^ 


183.4 ±22.1 


192.5 ±21.8 


185.4 ± 21.0 




-0.1 -0.05 0.05 0.1 0.15 0.2 -0.1 -0.05 0.05 0.1 0.15 0.2 

AE, GeV AE, GeV 



FIG. 3: The AE distributions for the B+ D+D^ (a) and the B~ D-D^ (b). The points with 
error bars show the data, the hatched histograms represent BB MC and the curves are projections 
from the 2D fits. 

given in Table [H where the hsted efficiencies include intermediate branching fractions. The 
projections of the 2D fit result is shown in Figs. [Hand [21 



Charge asymmetry in B^ D^D decay 



To calculate the charge asymmetry in the D^D^ decay channel, we repeat the 

fits separately for the D~^D^ and B^ D^D^ samples. The AE distributions for 

B~^ — > D^D^ and B~ —>■ D'D^ candidates are presented in Fig. [31 The fit results are given 
in Table [Til Using the results of the 2D fits, we calculate the charge asymmetry: 

N(D'D°) - N(D+D^) 
Acp = -4 ^ = 0.00 ± 0.08 1 

where the error is statistical only. 



Cross checks &; systematic uncertainties 



We calculate the B^ D^D^ branching fraction separately for all used D decay chan- 
nels, and the results are consistent with the average value. As an additional check, we apply 
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a similar procedure to a decay chain with similar final state: DfD^. We measure 

the branching fraction B{B^ D^ID^) = (0.95 ± 0.02)%, where the error is statistical only. 
This is consistent with the world average value (1.09 ±0.27)% [9|. The charge asymmetry in 
this final state is consistent with zero: (— 0.5±1.5)%. We also measure the charge asymmetry 
for the D^D^ background events and find a value consistent with 0: (—1.4 ± 1.3)%. 

The following sources of systematic errors are considered: tracking efficiency (6%), 
PID efficiency (2%), vr*^ reconstruction efficiency (6%), D branching fraction uncertainties 
(6%), signal yield determination (4%), luminosity determination (1.5%) and MC statistics 
(1%). The uncertainty in the tracking efficiency is estimated using partially reconstructed 
£)*+ — *• D^[Kg7r~^n^]7i^ decays. The kaon identification uncertainty is determined from 
D*~^ — >■ D'^[K~ 7r~^]7[^ decays. The error in signal yield determination is estimated by vary- 
ing the signal and backround shapes and fit range. We assume equal production rates for 
B^B~ and B^B^ pairs and do not include the uncertainty related to this assumption in 
the total systematic error. The overall systematic uncertainty is 10% for the branching 
fraction measurement. The systematic uncertainty is taken into account in the upper limit 
calculation: the signal efficiency is decreased by the systematic uncertainty. 

The asymmetry measurement contains the following systematic errors: tracking efficiency 
difference for vr^ (0.013), pion identification efficiency difference for vr^ (0.004) and signal 
yield determination (0.015). The total systematic uncertainty is 0.02. 



CONCLUSION 

In summary, we report improved measurements of B^ D^D^ branching fraction and 
the charge asymmetry: B{B+ D+D^) = (3.85 ± 0.31 ± 0.38) x 10"^ and Acp{B+ 
D+D^) = 0.00 ± 0.08 ± 0.02. We also set the upper limit for the B^ D^D^ decay 
branching fraction of B{B^ D^D^) < 0.42 x 10""^ at 90% CL. These results are consistent 
with our previous results [2] and supersede them. Our results are also consistent with BaBar 
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